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Abstract 



Feedback stabilization of an ensemble of non interacting half spins described by Bloch 

r \ , equations is considered. This system may be seen as a prototype for infinite dimensional 

systems with continuous spectrum. We propose an explicit feedback law that stabilizes 

asymptotically the system around a uniform state of spin +1/2 or -1/2. The proof of 

rS^ • the convergence is done locally around the equilibrium in the H^ topology. This local 

j^ I convergence is shown to be a weak asymptotic convergence for the H^ topology and thus 

a strong convergence for the C° topology. The proof relies on an adaptation of the LaSalle 

invariance principle to infinite dimensional systems. Numerical simulations illustrate the 

efficiency of these feedback laws, even for initial conditions far from the equilibrium 

T— I ■ 

>, 
'sj" ! Keywords: Nonlinear systems, Lyapunov stabilization, LaSalle invariance, ensemble con- 

ll^ ' trollability, infinite dimensional systems. 

^1: ; 1 Introduction 



1.1 Infinite dimensional systems with continuous spectra 



Most controllability results available for infinite dimensional bilinear systems are related to 
^ ■ systems with discrete spectra (see for instance, |2] for exact controllability results and [H |13] 

for approximate controllability results). As far as we know, very few controllability studies 
consider systems admitting a continuous part in their spectra. 

In [i ll] an approximate controllability result is given for a system with mixed discrete/continuous 
spectrum: the Schrodinger partial differential equation of a quantum particle in an N-dimensional 
decaying potential is shown to be approximately controllable (in infinite time) to the ground 
bounded state when the initial state is a linear superposition of bounded states. 

In [9l [To] a controllability notion, called ensemble controllability, is introduced and dis- 
cussed for quantum systems described by a family of ordinary differential equations (Bloch 
equations) depending continuously on a finite number of scalar parameters and with a finite 
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number of control inputs. Ensemble controllability means that it is possible to find open-loop 
controls that compensate for the dispersion in these scalar parameters: the goal is to simulta- 
neously steer a continuum of systems between states of interest with the same control input. 
Such continuous family of ordinary differential systems sharing the same control inputs can 
be seen as the prototype of infinite dimensional systems with purely continuous spectra. 

The article p.0| highlights the role of Lie algebras and non-commutativity in the design 
of a compensating control sequence and consequently in the characterization of ensemble 
controllability. In [3], this analysis is completed by functional analysis methods developed 
for infinite dimensional systems governed by partial differential equations (see, e.g., [8] for 
samples of these methods). Several mathematical answers are given, with discrimination 
between approximate and exact controllability, and finite time and infinite time controllability, 
for the Bloch equation. In particular, it is proved that a priori bounded L^-controls are not 
sufficient to achieve exact controllability, but unbounded controls (containing, for example 
sums of Dirac masses) allow to recover controllability. For example, it is proved in [3] that 
the Bloch equation is approximately controllable to 63, in H^, in finite time, with unbounded 
controls. The authors also propose explicit open loop (unbounded) controls for the local exact 
controllability to 63 in infinite time. 

The goal of this article is to investigate feedback stabilization of such specific infinite 
dimensional systems with continuous spectra. As in [11], the feedback design is based on a 
Lyapunov function closely related to the norm of the state space, a Banach space. 

1.2 The studied model 

We consider here an ensemble of non interacting half-spins in a static field {0,0, Bq) in M^, 
subject to a transverse radio frequency field {u{t),v{t),0) in M^ (the control input). The 
ensemble of half-spins is described by the magnetization vector M £ ]R.^ depending on time t 
but also on the Larmor frequency u = — 7-Bo (7 is the gyromagnetic ratio). It obeys to the 
Bloch equation: 

dM 

— {t, u) = {u{t)ei + i{t)e2 + wes) x M{t, oj), (1) 

where — cxd < uo^ < lo* < -|-oo, to G (w*, w*), (ei, 62, 63) is the canonical basis of M'^, x denotes 
the vector product on M^. The equation ([!]) is an infinite dimensional bilinear control system 
in which 

• the state is the function M, with, for every u £ (a;^,,tj*), M{t,io) G S^, the unit sphere 
of R^ 

• the two control inputs u and v are real valued. 

It must be stressed that u{t) and v{t) are common controls for all the members of the ensemble, 
and they cannot depend on w. In coordinates M = [x, y, z), the Bloch equation may be written 



(2) 



where x stands for dx/dt. 

Now, let us precise the definition of a solution associated to Dirac controls. When u,v £ 
-^Lc(^)' then, for every initial condition Mq £ L^((a;*,a;*),M^), the equation ([1]) has a weak 



X = 


= —ujy + vz. 


y = 


= UJX — uz. 


z = 


= —vx + uy. 



solution in the usual sense: M G C*^([0, +oo),L^((a;:^,a;*),M'^)). Let us give the definition 
of a solution associated to a Dirac control. In the sequel 5a{t) stands for the Dirac function 
5{t — a). If u = a5a{t) + u^ and v = V'^ where u^^v^ S Lj^J^)^ a > and a G (0, +00), then 
the solution is the classical solution on [0, a) and (a, +00), it is discontinuous at the point a, 
with the explicit discontinuity 

/I \ 

M(a+,a;)= cos(a) - sin(Q) \M{a~,oj). 
\ sin(Q) cos(q) / 

This definition corresponds to the limit, as e — t- 0^ of solutions associated to n(t) = {a / e)\\a^a+e]{t) ■ 
Formally, the spectrum of the operator A defined by 

{AM){oj) := wes x M{uj) 

is — i(w*,a;*) U i(a;*,a;*): for every wjj G {lo^^lo*), the eigenvector associated to itiwu is 
{1,^1,01)^ 5 i^Auj). Thus the Bloch equation ([T|) is a prototype of infinite dimensional system 
with continuous spectrum. 

1.3 Outline 

The goal of this article it to propose a first answer to the following question. 

Local Stabilization Problem. Define an explicit control law {u,v) = {u{t,M),v{t,M)) 
and a neighborhood U of —63 (in some space of functions (a;*,a;*) — )• S^ to be determined) 
such that, given any initial condition M" S U, the solution of the closed loop system is defined 
for every t E [0, +00), is unique and converges to —63, when t — )• +00, uniformly with respect 

to w G (a;*,a;*). 

Section [2] is devoted to control design and closed- loop simulations: the feedback law is 
the sum of a Dirac comb and a time-periodic feedback law based on a Lyapunov function; 
Proposition [1] proved in Appendix A guarantees that the closed-loop initial value problem is 
always well defined; simulations illustrate the convergence rates observed for an initial state 
formed by a quarter of the equator on the Bloch sphere. In section [3] we state and prove the 
main convergence result. Theorem [1) the closed- loop convergence towards the constant profile 
M{u)) = —63 is shown to be local and weak for the H^ topology on M. The obstruction 
to global stabilization is also discussed: it is based on an explicit description of the Lasalle 
invariant set. Some concluding remarks are gathered in section U) 

2 Lyapunov H^ approach 

2.1 The impulse-train control 

It is proved in [3] that controls containing sums of Dirac masses are crucial to achieve the 
controllability of the Bloch equation. In view of the controls used in this reference, it is 
natural to consider a control with the following "impulse-train" structure 

+00 
u = u + ^Tr 6iir{t), v = {-lf(^')v (3) 

k=l 



for some period T > 0, which is fixed in all the article {E (7) denotes the integer part of the 
real number 7). The new controls u and v belong to Lj^^{W). Then, after each impulse that 
is applied at time t = kT, x remains unchanged, but y and z are moved to their opposites, 
that is 

{x, y, z){kT+) = (x, -y, -z){kT'). 

The resulting state diffeomorphism 

{x,y,z)^{X = x,y = -y,Z = -z) (4) 

transforms ([2|) into 

X{t) = u;y{t)-v{t)Z{t), 
y{t) = -ujX{t)-uZ{t), 

z{t) = v{t)x{t) + u{t)y{t). 

Let q = (—1) yTj ^ Considering the following change of variables 

{x,y,z){t,u) ^ {^x{t,uj),Cit)y{t,u}),Cit)z{t,uj] 

one gets the following dynamics 

X = —<^u;y + vz, 

ij = <;ijjx — uz, (5) 

z = —vx + uy, 

with the new control {u,v) as in ([3|). It is as if, between [kT, (fc + l)T] and [{k + l)T, {k + 2)T], 
one is changing the sign of cj, but the solution, after the identification (j4]), remains continuous 
in t (but not differentiable in t at the instants t = kT, k G'N). In other words, the application 
of the impulses at i = kT changes the sense of rotation of the null input solution. One 
would expect that this impulse-train control is reducing the average dispersion of the solution. 
Roughly speaking, the dispersion observed for the open-loop system ([2]) with (u, v) as input 
is strongly reduced and almost canceled for the open-loop system ([5]) with (u, v) as input. 

2.2 Heuristics of the Lyapunov-like control 

Now let Z{t,uj) and Q{t) defined by 

Z = x + iy, n = V — iu 
where x,y,z refer to the transformed dynamics ([5]). Then one may write ([5]) in the form 

iq{t)ujZ{t,uj) + Q,(t)z{t,uj), 




-^ 



n{t)z{t,uj) 



where 3f?(^) (resp. ^) denotes the real part (resp. the complex conjugate number) of a complex 
number ^ G C. The following transformation 

Z{t,uj) = Z{t,uj)e-'^-fo<^)'i^ 



converts the system into the driftless form 



Z{t,co) = n{t)z{t,Lo)e-^'^^o^, 



z{t,uj) 



-R 



n{t)Z{t,u)e-'^-foi 



(6) 



where, for notation simpUcity, one lets Z{t,uj) stand for Z{t,uj), and one lets L <^ stand for 

/o <i{r)dT. 

For the moment one shall assume that the input il.{t) will be chosen in such a way that 
the solution {Z{t,uj),z{t,u})) of ([6]) does exist, it is unique and it is regular enough in a way 
that one may consider that the derivatives Z'{t,uj) = ■gj(i,w) and z'{t,Lo) = Q^{t,uj) exist 
almost everywhere and they are solutions of the differential equation that is obtained by 
differentiation of ([6]) with respect to u, namely 



z' = n 



z! 



-^ 



z' -i{.fn'^)Z 



'k 



n 



^]z 



0^}, 



-}' 



where Z' stands for ^Z' , and z' stands for ^z'. 



- dt^ ■ 
Now consider the following Lyapunov-like functional: 

£ = 1 f {\zf + (/)2 + 2G(z + l)}du; 



(7) 



(8) 



where G is a positive real number and Z{t,uj), Z'{t,uj) and z'(t,oj) refer to the solutions 
respectively of ([6]) and ([7|). One may write 



jC{t) = ^[ I [Z'Z' + z'z' + Gi} duj 



(9) 



and so, taking into account ([6]) and ([7]), the fact that rt{t) does not depend on w, one gets 

d 



dt 



C{t) = R[n{t)H{t)] 



(10) 



where 



Hit) :- 



J {iijA [Zz'-Z'z)-Gz]i 
w, \o / 



'k'o^duj. 



Hence one may take Vt{t) = —KpH(t), where Kp is a positive real number, obtaining 
n{t) = Kp Hil [A {Zzl - Z'z) + Gz}e*'^^o ?da;. 



:iii 



oj, \o 



It follows that 



^<'> = -]^l"«i^- 



Consider the system ([6]) in closed loop with the control law (jlip . thereby called by closed 
loop system. The state space of this system is /f-^((a;^,, a;*),M'^), which is the set of functions 
/ S L'^{uj^,uj*) such that the distributional derivative /' belongs to L^(w*,w*). This space, 
equipped with the norm 



1/2 



H^ 



\f'{u;)\' + \f{u)\'dco 



is a Banach space. In other words, closed-loop system (l6])- (llip may be considered to be a 
differential equation of the form 

M{t) = F{t,M{t)) 

M(0) = M^ e H\{uj^,uj*),S'^) 

where F{t, M) is a continuous map 

F : M X //i((a;„w*),M3) -^ H^{{uj^,u*),^^). 

Moreover, F is locally integrable (Lj^^) and periodic in t and locally Lipschitz in M. Using 
the same ideas as in the proof of the Cauchy-Lipschitz (Picard-Lindelof ) theorem, we get local 
(in time) solutions in C^{[0,T], H^). From the contruction of the feedback law, finite time 
blow up in H^ is impossible, thus solutions are global in time. Precisely, we have the following 
result, whose proof is detailed in Appendix A. 

Proposition 1 For every initial condition M^ E if"^((a;*, u;*),§^), the closed loop system ^, 
[Tl\) has a unique solution M £ C^^ ([0,oo),i^^ ((a;*,a;*),S2)) such that M(0) = M°. 

In this statement, the space Cp^ [[0,oo), H^ ((a;*,CiJ*),S^)) is made of functions 
M E C°([0,oo),Fi ((a;,,a;*),S2)) such that M E C^{{kT,{k + 1)T),H^ {{uj^,uj*),§'^)), for 
every /c E N; their derivative dM/dt is possibly discontinuous at i = kT, but it has finite 
limits in H^ {{uj^,uj*),§,'^) when t -^ {kT)+ and t -^ {kTy . 

2.3 Closed-loop simulations 

We assume here cj^, = 0, w* = 1 and we solve numerically the T-periodic system ([5]) with the T- 
periodic feedback law (jlip where Z = x + iy and Vt = v — iu. The parameters are T = 2'K/{ijj* — 
w^), Kp = l,G = rV20. The simulation is for t E [0,r/], T/ = 20T. The w-profile [w*,w*] 3 
uj I—)- {x{t,uj),y{t,uj), z{t,uj)) is discretized {0, . . . , N} 3 k ^ {xk{t),yk{t), Zk{t)) with a regular 
mesh of step ejsi = '^ J^'^* and A'^ = 100: {xk{t),yk{t), Zk{t)) is then an approximation of 
{x{t,k€N),y{t,k€]sr),z{t,kei\f)). We have checked that the closed-loop simulations are almost 
identical for A^ = 100 and A^ = 200. In the feedback law (jlip . the integral versus u is computed 
assuming that {x,y,z) and (x',y',z') are constant over ]{k — ^)eN,{k + ^)e7v[, their values 
being {xk,yk,Zk) and ( '"'*2^^'''\ ^'°^2~eT'^ ' '"'^2rw'"0 ' "^^^ obtained differential system is of 
dimension 3(A^ +1)- It is integrated via an explicit Euler scheme with a step size h = T/800. 
We have tested that h = T/1600 yields to almost the same numerical solution at t = Ty = 20T. 
After each time-step the new values of {xk,yk, Zk) are normalized to remain in S^. 

Figures [1] and [2] summarize the main convergence issues when the initial w-profiles of 
{x,y,z) E §2 are z{0,uj) = 0, x(0,a;) = cos (f w -F f ) and y(0,cj) = sin(|a; + f). The 



convergence speed is rapid at the beginning and tends to decrease at the end. We start 
with £(0) ss 3.24 and sup^^j,^^^,^,] ||M(0,cj) + esU w 1.41. We get C{20T) ^ 0.38 and 
sup^gf^^ ^*i ||M(20r, w) + esll w 0.27. The control problem is quite hard due to the fact 
that one has a continuous spectrum, that is, an infinite ensemble of systems with a com- 
mon control input J7(t). Hence, as time increases, the control must fight against the dis- 
persion of the solutions M(t,uj) for different values of u. Simulations (not presented here) 
on much longer times until lOOOT and with the same initial conditions and parameters al- 
ways yield to smaller final value for the Lyapunov function: we get £(1000T) ~ 0.02 and 
sup^gr^, 1^*1 ||-^(1000T, Li;) + esll w 0.04. This is a strong indication of asymptotic converge 
of the profile M(t,uj) toward —63, even if the convergence speed seems to be very slow. This 
numerically observed convergence is confirmed by Theorem [T] here below. 
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Figure 1: Lyapunov function C{t) defined by ([8]) and the closed- loop control il.{t) defined 
by (HH) . 
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Figure 2: Initial [t = 0) and final (t = Tf) tj-profiles for x, y and z solutions of the closed-loop 
system ([5]) with the T-periodic feedback (jlip . 



3 Main Result 

3.1 Local stabilization 

The main result of this paper shows that the control law (jlip is a solution of the local stabi- 
lization problem stated at the end of the introduction. 

Theorem 1 Consider system ([6]) with the feedback law ([TT]) . There exists 6' > such that, 
for every Mq G //-^((o;*, a;*),§^) with \\Mq -\- eaU/^i < 6', M{t,ijj) converges weakly in H^ to 
—63 when t —7- -|-oo. In particular, as the injection of H^ in C^ is compact, M{t,uj) converges 
to —63 when t — )• -|-oo uniformly with respect to uj £ (a;*,a;*) (convergence in the sup norm of 

The proof of Theorem[T]relies on an adaptation of the LaSalle invariance principle to infinite 
dimensional systems. This principle is a powerful tool to prove the asymptotic stability of 
an equilibrium of a dynamical system in finite dimension: one just needs to check that the 
invariant set coincides with the target [12]. For infinite dimensional systems, the use of the 



LaSalle invariance principle is more delicate (because closed and bounded subsets are not 
necessarily compact). Roughly speaking, there are 2 adaptation strategies: 

• either one accepts a weaker result: approximate stabilization dKH] or weak stabilization 
[1]; in this case, one may only need to ensure that the invariant set is reduced to the 
target, 

• or one wants strong stabilization; in this case, one needs an additional compactness 
property for the trajectories of the closed loop system [6], or a strict Lyapunov function 

m- 

This article concerns the first strategy: we prove weak stabilization. 

The first step of our proof consists in checking that, locally, the invariant set is {—63}. 

Proposition 2 There exists 6 > such that, for every M° G if^((a;*,cj*),S^) with ||M° + 
^sll^i < ^! ^^6 map 1 1—)- C{t) is constant on [0, +00) if and only if M^ = —63. 

Proof: Let us assume that £(i) is constant. Then, r2(t) = 0, Z(t,u}) = Zq{oj) and z(t,a;) = 
zq{lo). We deduce from (llip that 

it{Z'oZo - Zqz'o) - GZo\e''^^du; = 0,Vt G [0,r]. 
Considering the power series expansion versus t of the left hand side, we get 



P'iZ'f^zo - ZM - GPZq duo = 0,VP G CM. (12) 

Polynomials are dense in H^{uj^,uj*), thus, the previous inequality holds for every P G 
H^{lo^,uj*). In particular, with P{io) = Zq{lo), we get 

Zl(^Z'oZo - Zoz'o) - GZ;;Zo\ du = 0. 
Summing this equality with the following left hand side, we deduce that 



Juj 



Z^\' + G\Zo 



doj 



{I + zo)\Z'qV - Z'^Zoz'q 
2 

L2 



^ II 1 + ^o||l°° 11-^011x2 + 11^0 11^2 



dbJ 

7 II lU' II 

^o\\L°° IFoIIl2. 



Thanks to the continuity of the embedding H^{uj^,uj*) C L°°(a;*,a;*), there exists a constant 
C > (independent of M) such that 



1/2 



^oIIl- ^C{ {\Z'o\' + G\Zo\')doo 



Therefore, we get 



/:; iz^p+G|Zoi 



70 P dio. 



^ [\\1 + zo\\l^+C\\z'Jl-J^ [\Z'^\^ + G\Zo 

There exists 6 > such that, for every Mq G H^{{uj^,,co*),S'^) with \\Mq +es\\i{i < 6, we have 

||1 + zo\\l°° +C\\z'f)\\i2 < 1. 

If C is constant along the trajectory associated to such an initial condition Mq, then, the 
previous argument shows that Zq = 0, thus Mq = —63. D 

For the proof of Theorem [H we need the continuity with respect to initial conditions, of 
the solutions of the closed loop system ([6]), (jlip . for the ff 2 (a;,,,a;*)-topology. This space is 
defined by interpolation between L^(a;^,a;*) and H^{io^,UJ*) and we have a compact injection 
H^{lo^,uj*) —7- H^{co^,uj*). We also use the space H~^{uj^,ui*), which is the dual space of 
H2(^uj^:,u)*) for the L^-scalar product. First, let us recall the following Lemma. 

Lemma 1 There exists ci > such that, for every if G H2(^u}^,u}*) and for every a G [0,T], 
the map co 1— t- (/7(a;)e*"'^ belongs to i?2(u;^,,w*) and satisfies 

LO)e 1 < Cl 03 1 . 

Proof: We have 

||(^He^"'^||i2 = ||¥p||i2,V(^ G L2(a;„w*), 

and, for every 99 G ff"^(u;*,cj*), 

||(^(a;)e^°'^||^i = / y {iS) + ia^{i^)\^ + \^{i^)\^du: 
^ / 2|v3'(a;)|2 + (2a2 + l)\ip{uj)\^ du) 

thus we get the conclusion with, for example, ci := (2T^ + 2)^''^ by interpolation. D 

Proposition 3 There exists 5' > such that, for every (M°)„gN G i:/'H(^*)'^*))S^)^7 ^Sd £ 
H^{{ui^,(jj*),'B^) satisfying 

• ||M0 + 6311^1 <<5',VnGN, 

• M^ ^ M^ weakly in H^ when n — )■ +00, 

• M^ —7- M^ strongly in H2 when n — t- +cxd, 

i/ie solutions Mn{t,u}), M^{t,uj) of the closed loop system associated to these initial conditions 
satisfy the following convergences, when n — )• +00, for every t G [0, +00), 

Mn{t) — ^ Moo{t) strongly in H2 and $7„(i) — )• S7oo(i)- 

10 



Proof: First, let us emphasize that y^C{M) and ||M + e3||/^i are equivalent norms on a small 
enough //^((ci;^,,cj*),§^)-neighborhood of —63: there exists 7/, c*,c* > such that, for every 
M E H'^{{uj^,uj*),S'^) with ||M + e3||^i < i], we have 



c,y^C{M) < \\M + 631^1 < c*^/C{M). 



(13) 



Now, let 6' := mm{SC:^/c* ,t]}, where 6 is as in Proposition [2l Thanks to the monotonicity of 
C, we have, for every t £ [0, +00), 



We have 



\Mn{t) +6^11^1 <C*y /C{Mn{ t)) 

<^||M0 + e3ki<^<<^. 



\Mnit)-M^ 



+ f^\\F{s,Mn{s))-F{s,M^{i 



1 < ||M0-M0 II 1 
I 1 ds. 



Let us prove the existence of C > such that, for every M,M G H^{u^,,lj*) satisfying 
||M + e3||j|/i < 6, we have 

\\F{s,M)-F{s,M)\\^i <C||M-M||^i,VsGM. 

Then, we will conclude the proof thanks to the Gronwall Lemma. Let us work, for example, 
on the third component of F: 

Fs{t,M) = -Re \n{t)Ze-''^-^o( 

where Q is defined by (llip . We have 

||F3(t,M)-F3(t,M)|| 



1 



< |l^(t)- ^(t)|||Ze — /oqi^i 
+|f2(i)|||(Z-Z)e-*-/oC| 

< \nit)-m\ci\\z\\^^+Kci\\z-z\\^. 



HI. 



1 



where ci is as in the previous Lemma and K = K{5). It is sufficient to prove the existence of 
a constant C > such that, for every M,M e H^{{uj^,uj*),E>'^) satisfying ||M + e3||j:/i, ||M + 
63||_ffi < ^, we have 

\n{t) - ^{t)\ < C\\M -M\\ 1 , Vt G [0, +00). 

H^ 
Let us prove it only on one of the terms that compose il (the other terms may be treated as 
well): 

Jf^ (z'z-Z'z\e''^-foiduj 
< jfiz'-Z'^e'^^^oCdu 
+ Jf Z'(z-z]e''^Io(doj 



< 


\\Z'-Z'\\ illze^'^Jo'^ll 1 




+\\z'\\ 111(^-5)6*^/0 "^11 1 


< 


ciK||M-M||^i. D 



11 



Proof of Theorem [1} Let 6' be as in the previous proof. Let M^ G H^{{uj^,uj*),S'^) be such 
that ||M° + e3||j:/i < 6' and M G C^{[0,+oo),H^{{uj^,uj*),8'^)) be the solution of the closed 
loop system such that M(0) = M^. 

First step: Let us prove that J7(t) — t- when t — t- +cxd. Thanks to the choice of the feedback 
law, M{t) is bounded in H^, uniformly with respect to t S [0, +oo). Computing explicitly 
-^(t), we see that -^(t) is bounded in C uniformly with respect to t G [0, +oo) — NT. Thus, 
Q is uniformly continuous on [0, +oo). Since fi G T^(0, +oo), it has to satisfy $7(f) — )■ when 
t — )• +00 (Barbalat's lemma). 

Second step: Let us prove that —63 is the only possible weak H^ limit. Let M^ be a weak 
H^ limit of the trajectory starting from M^. There exists a sequence (i„)ngN of [0, +00) such 
that tn — >• +00, 

M{tn) -^ M^ weakly in H^ when n — > +00, 
M{tn) — )■ M^ strongly in if ^/^ when n — )- +00. 

Working as in the previous proof, one may prove that 

l|Af(t„) + e3||^i <5,VnGN. (14) 

There exists t^o £ [0, T) such that t„ mod T — )• t^o- Let M^{t,uj) be the solution of the 
closed loop system associated to the initial condition Moo(icxD) = ^^- Let us prove that C is 
constant along this trajectory, by proving that the associated control Vt^ vanishes. In order 
to simplify, we assume that ioo = (otherwise, consider an additional shift). For every t > 0, 
M{tn + t) — )• Moo(t) strongly in H^'^ when n — )• +00, thanks to the previous proposition. 
This allows to pass to the limit in the feedback law: il(t„ + t) — )• ^ooit) when n — )• +00, for 
every t > 0. Thanks to the first step, we get Q^o = 0. 

In order to apply Proposition [21 we only need to check that ||M^ + esU/^i < 6, which is a 
consequence of ([T4|) . D 



3.2 Obstructions to global stabilization 

Now, let us explain why these feedback laws may not provide global stabilization in H^{(uj^,uj*),S'^) 
The first obstuction is a topological one: the space ii^((a;^,,a;*),S^) cannot be continuously 
deformed to one point (because S^ is not), thus global stabilization in this space is impossible. 
Actually, for our explicit feedback laws, it is easy to see that M^ = +63 is an invariant 
solution. It is interesting to know whether it is the only one (i.e. if one may expect the 
stabilization of any initial condition M^ 7^ 63). The answer is no, as emphasized in the 
following proposition. 

Proposition 4 For every U}^,u}* G M such that w* < uj* , there exists an infinite number of 
non trivial functions in the LaSalle invariant set. 

The proof is detailed in Appendix[Bl Actually, all the invariant solutions may be computed 
explicitly. 
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4 Conclusion 

We have investigated here the stabihzation of an infinite dimensional system admitting a 
continuous spectrum. We have designed a Lyapunov based feedback. Closed-loop simulations 
illustrate the asymptotic convergence towards the goal steady-state. We have provided a local 
and weak convergence result for the H^ topology. Simulations indicate that the domain of 
attraction is far from being local and thus we can expect a large attraction domain for this 
feedback law. However, the stabilization is not global because there exists non trivial invariant 
solutions. 

Few problems are still open concerning this problem. Are the invariant solutions unstable? 
Does the local stabilization hold for the strong i7^-topology (not only the weak one)? Is it 
possible to get semi-global stabilization? What is the value of convergence rates? Is it possible 
to produce arbitrarily fast stabilization? 

More generally, this feedback and convergence analysis opens the way to asymptotic stabi- 
lization of neutrally stable systems of infinite dimension with continuous spectra. For example, 
it will be interesting to see if the following system (ID Maxwell-Lorentz model for the prop- 
agation of an electro-magnetic wave in a non-homogeneous dispersive material) can also be 
stabilized to zero: 

^=p^{x){E-P), XG(0,1) 
E{0,t)=u{t), E{l,t)=v{t) 

with two controls n and v. When p{x) is a smooth strictly increasing positive function, the 
above system admits as continuous spectrum zizi]p{0) , p(l)[. 
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A Proof of Proposition [1] 

Let M° G i?i((a;*,a;*),S2) and ii > be such that 

/2>max|||M°||j:^i,V^2£(0)-Fw* -a;*| (15) 

Let Ci , C2 > be such that 



\\fe-'^'\\m<C4f\\HU 

yf eH^{uj^,uj*),yte [o,T], 

\\F{t,M^)-F{t,M2)\\m <C2\m-M2\\m, 
VMi,M2 G Sfi[Fi((w„w*),M3)],Vt e [0,T], 



(16) 



(17) 



where i?j:j[X] denote the closed ball centered at with radius R, of the space X. Let T* = 
T*{R) > be small enough so that 

\\M°\\hi + T*CiKp{G + 2T)R^ < R and T*C2 < 1. (18) 

Let us consider the map 0, defined on the space 

E := Bn[C\[0,T*],H\iu,,u*),R^))] 
by 

e(M)(i,a;) :=M°(w)+ / F{s,M{s,uj))ds 



JO 

for every {t,uj) G [0,T*] x (a;*,w*). 

First step: Let us prove that Q takes values in E. Let M €z E. It is clear that G(M) is 
continuous in time with values in H^{{uJ:^,uj*),'R^). For t £ [0,T*], we have 

||G(M)(t)||^i < ||MO||^i + / \\Fis,MismHids. 

Jo 
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By definition, we have 



\\F{s,M{s))f^, 
= \\n{s)z{s)e-'^fo^\\j^^ + P[0(s)Z(s)e-*^/o ^]||^, 

< \n{srcf{\\z{s)\\l, + \\z{s)\\l,) 
= cf\n{snM{s)\\j,,. 

Moreover, the Cauchy-Schwarz inequahty gives 

Ms)\<Kp{G + 2)\\M{s)\\j,„ 
thus, thanks to (|18p . we have 

||0(M)||ioo((o,r*),//i) 
< \\M'^\\Hi+T*CiKp{G + 2)R^ <R. 

Second step: Let is prove that @ is a contraction. For Mi,M2 G E and t G [0,T*], using 
dni), we get 

||e(Mi)(t)-e(M2)(t)|bi 

< tC2||Mi - M2||ioo((0,T*),Hi)) 

thus is a contraction, thanks to (llSp . 

Third step: Let us prove the existence and uniqueness of strong solutions, defined on 
[0, +00). Thanks to the Banach fixed point theorem, the map O has a unique fixed point. 

We have proved that, for every R > 0, there exists T* = T*{R) > such that, for every 
M" G Bji[H^{{uj^,,uj*),S'^)], there exists a unique weak solution 

M G C\[0,T*],H\{uj^,u*),M.^)) 

in the sense 

M{t,uj) = M'^iij) + J^F{s,M{s,u}))ds, 
mH^{io^,io*),yt€ [0,T*]. 

From this equahty, we deduce that 

and 

^{t,u) = F{t,M{t,u;)) 
in H^{{uj^,uj*),R^)yt£ [0,T*] -NT. 

Since H^ C C^, we also have 

^{t,io) = F{t,M{t,Lo)) 

VtG [0,r*]-NT,Va;G K,a;*). 

This has 2 consequences: 

• M(t, .) takes values in S^ for every t G [0,T*], indeed, M^ does and the following 
computation is licit for every t G [0, T*] — NT, 

j^\\M{t,u)f = 2{M{t,u),F{t,M{t,u))) = 0, 
15 



• the computations ([9]), (jlOp are licit, thus C{t) is not increasing. 
Therefore, we have 

mT*)\\i, 

= / |z'(r)|2 + z'(r)2 + \z{T)\^ + z{Tfduj 

< 2£(T) + a;* - w* 

< 2£(0) + w* - w* < i?2 

thanks to (|15p . Thus, we can apply the previous result with M^ replaced by M{T*): it provides 
a solution on [0, 2T*]. Iterating this again, we get a solution defined for every t G [0, +oo). D 

B Proof of Proposition |4] 

B.O.l Heuristic 

In this section, we perform an heuristic to prove that the LaSalle invariant set contains non 
trivial solutions. In order to simplify the notations, we take G = 1. 

Let M = (x, y, z) be in the LaSalle invariant set. In view of (|12p and after integrating by 
parts, we have 

xz" — zx" = X on [bj^^uj*), 

yz" - zy" = y on (l^*, cj*), ,^g- 

x^ + y'^ -\- z"^ = 1 on {(jj^^uj*), ^ ' 

xz' = x' z and yz' = y'z at u}^,,u}*. 

First step: We proceed to eliminations in order to get an ordinary differential equation 
involving only z. Differentiating the third equality of (J19p . we get 



xx' + yy' + zz' = on (a;*,a;*), (20) 

xx" + {x'f + yy" + {y'f + zz" + {z'f = on iu;,,u*). (21) 
Multiplying this equality by z and using the 3 first equalities of (J19p . we get 

z" + z^-l + z[{x'f + (y'f + {z'f] = on (^„l^*). (22) 



Derivating this equality, using again the 2 first equalities of (J19p together with (|22p and (|20 
we get 

zz'" + 3z^z' + z'{l-z") =0on (a;*,w*). 



At the points cj=k and to*, thanks to (j20p and the 4th equality of (J19p . we have 

= z(xx' + yy' + zz') = z'ix"^ + y^ + z^), 

thus 

z' = at w* and w*. 

We deduce from the 4th equality of (|19p that 

zx' = zy' = at cj=K and uj*. 
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Thus, the relation (j22p provides 

z" ^ z — 1 = at cj^, and w*. 

Therefore, if M = (x, y, z) is in the invariant set, then z solves the following boundary value 
problem 

zz'" + -iz^z' + z\\ - z") = on {uj^,oj*), 

z' = z" + z^ -1 = at uj^ and w*, ^ ' 



and (x, y) are solutions of the first order system 

xx' + yy' = —zz' on (w*, w*), 

[x'f + {y'f = i^4^ - {z'f on (a;., a;*), (24) 

zx' = zy' = at a;*,a;*. 

Second step: Let us solve the equation (I23p . We introduce the function F := z" + 82;^. The 
first equality of (|23p allows to prove that 



zF' = {F -l)z on (w*,a;*). 
Thus, there exists C E M such that [F — 1) = Cz, i.e. 

z" = -3z'^ + Cz + lon{u}^,u}*). (25) 

Thanks to this equation, we deduce from the second equality of (p3|) that 

z{C — 2z) at w* and w*. (26) 

Multiplying ([25]) by z' and integrating over {uj^:,uj), we get 

-(z'f = -z^ H z'^ +Z + est. 

2^ ^ 2 

The left hand side vanishes at the boundary, and the right hand side is equal to z at the 
boundary thanks to ([26|) . thus 

Z{UJ^:) = Z{U)*). 

The conclusion of this second step is the existence of a constant C G M such that z solves one 
of the following systems 

fS )l ^" ^ -Sz^ + Cz+l on (w*,w*), 
'^ ^' y z = z' = at uj^,uj*, 

z" = -3z^ + Cz + 1 on (w*,a;*), 
z — C/2 = z' = at w*, w*. 

Third step: We prove that (Si) has admissible solutions for arbitrarily small intervals 
{uj^:,uj*). Multiplying the first equality of (Si) by z' and integrating over (w*,w), we get 

{z'f = -2z^ + Cz^ + 2zon{uj„uj*). (27) 
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The function t H> -2*^ + Ct'^ + 2t vanishes at t = 0, t = (3c ■= (C + VC^ + 16)/4 (simple 
roots) and is positive on {0,/3c) thus one may define 

X 



ac := G ^ j . 

Then G G C°^(0,/3c), G is increasing from to ac and it has an infinite derivative at x = 
and X = j3c- Thus, G^^ : [0, ac] -^ [0, /3c] belongs to C'^{0,ac), it is increasing from to /3c 
and its derivative vanishes at and ac- Then the function z : (0, 2ac) ~^ I^i symmetric with 
respect to ac and such that 

z{u}) := G-^{uj),yuj £ [0,ac] (28) 

is a solution of (Si) with w* = and oj* = 2ac- Let us emphasize that this solution is 
admissible, when C < because it takes values in [0, /3c] which in included in [0, 1). Notice 
that 

ac > 0, 

thus, we have built admissible solutions of (Si) for arbitrarily small intervals (a;*,a;*) = 
(0,2ac). 

Fourth step: We prove that, for any solution of (Si) with C < 0, the system (I24p has 
solutions. First, notice that, thanks to (Si) and ([27|) . we have 



1 - z2 _ z" 



'\2 



2z-^ - Cz^ - C. 



z 
Eliminating x' in the two first equalities of (pi]) , we get 

a{y'f + 6y' + c = 

where 

a := x2 + y2 = 1 _ ^2^ 

h := 2zz'y, 

c := {zz'y - x'^{2z^ - Gz^ - C). 

Thanks to (I27p one may prove that the discriminant is 

A:=52-4ac = -4Cx2. 
Therefore, 



By symmetry, we also have 



y' = 




x' = 


-^' ,.^\v 



l-Z^ l-z2- 
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In order to ensure xx' + yy' + zz' = 0, the signs it need to be opposite. Therefore (x, y) are 
solutions, for instance, of the following linear system 



1-2 



rX + 






y = -t^y - f^x, (29) 

x(0) = xo, 

, y(0) = yo, 

where xq and yo a-re I'eal numbers such that Xg + yg = 1. 

Conclusion: For every C < 0, we have built candidates of invariant solutions of the closed 
loop system associated to the interval (a;^,,a;*) = (0,20^7), which is arbitrarily small when 
C — )• —00. In order to conclude, one just needs to check that this candidate indeed solves 
(J19p . which will be done rigourously in the next subsection. 



Remark 1 Let us emphasize that the same phenom,ena happens if we put different gains in 
the Lyapunov functions: for every Gi,G2 ^ 0, the feedback laws associated to the control 
Lyapunov function 

Cit) := / \Z'f + {z'f + Giz + G2\Zf du 



generate a non trivial LaSalle invariant set. 



B.0.2 Rigorous proof 

Lemma 2 There exists a continuous function 



(-00,0) -> (0, +00) 
G ^ Lc 



such that, 



for every G G (—00, 0), there exists a function z £ C^([0, Lc], [0, 1)) such that z > on 
{0,Lc) and 

' z" = -3^2 + Gz + l on (0, Lc), 
z = z' = at and Lq, 



Lc —7- when G 



-00. 



Proof: Consider Lc '■= 2ac and z : [0, Lc] — )• M, symetric with respect to ac and defined by 

dSHD- □ 

Proof of Proposition m Let -00 < oj^ < 00* < +00. Let {N,G) G N* x M;^ be such that 
w* — a;* = NLc (the existence is ensured by the intermediate values theorem and there exists 
an infinite number of such couples). Let z : [0, Lc] — )• [0,1) be as in the previous Lemma. 
Let Xo,yo S M be such that Xg + yg = 1, and (x,y) be the solution of (f29|) : this solution is 
well defined on the whole interval [0, Lc] because the system is linear and its coefficients are 
continuous. Now, let us check that (fT9|) holds with (a;*,a;*) replaced by (0,Lc). 
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First step: We check that the 3rd equahty of (fT9|) holds. Thanks to (|29|) . the quantity 

N := x^ + y'^ + z^ solves 

n 7 7 

-(l-iV) = ^— -^(l-iV)on(0,Lc) 

and (1 - iV)(0) = 0, thus iV = 1 on (0, Lq). 

Second step: We check that the fourth equality of (fTOl) holds. Since z = z' = at and Lc, 
we also have xz' = x'z and yz' = y'z at and Lq- 



Third step: We check that the 2 first equations of (J19p hold. The computations are similar 
to the ones of the Heuristic but now, we know that the functions considered are smooth (they 
are explicit), so these computations are licit. Starting from (I29p . using the result of the first 
step, the equality (|27p . and the 2nd order equation solved by z, we get 

{x'f + {y'f + {z'f = 2z-C= iz.i^. (30) 

Thus, we have 

z" + z''-l + z[{x'f + {y'f + {z'f] = 0. (31) 

The result of the first step justifies 

{x' + y^ + z')z" + [{x'f + {y'f + {z'f]z - {x^ + y^) = 0, 

that may be written 

x{xz" -x)+ y{yz" - y) + z'^z" + [{x'f + {y'f + {z'f]z = 0. (32) 

Now, differentiating 2 times the identity x^ + y^ + z^ = 1 and multiplying the resulting equality 
by z, we get 

xx"z + yy"z + z^z" + [{x'f + {y'f + {z'f]z = (33) 



Thus, dMD-dSD gives 

x{xz" — X — x" z) + y{yz" — y — y" z) = 0. (34) 



Derivating (I3ip and multiplying the resulting equality by z, we get 

zz'" + 2z^z' + z'{\ -z^ - z") + 2z\x'x" + y'y" + z'z") = 0. (35) 

We deduce from the second order equation solved by z that 

zz'" + 3z^z' + z'{l - z") = 0. (36) 

Indeed, we have 



zz 



ji' 



z(-6zz' + Cz') 



Sz'^z' + z'{l - z") = z'[3z^ + l + 3z^ -Cz-l]= zz'{6z - C). 
Thanks to the identity xx' + yy' + zz' = 0, the equation (l36|l may be written 

zz'" + 2z^z' + z'{\ -z^ - z") + 2z[z"{xx' + yy' + zz') + zz'\ = 
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or 

ZZ'" + IZ'^Z' + Z'(l - z2 _ z") 

+2z[x'{z"x -x) + y'{z"y - y) + z"z'z] = ^ ^ 

Finally, ([STD-dSSD gives 

x'{xz" -X- x"z) + y'{yz" -y- y" z) = 0. (38) 

Now, (j34p and (|38p give the conclusion, because 



X y A V |C| , 2 , 2^ 
x' y' J 1 - z2 



det( :, :, )=^^^(2;^ + y') = \/|C|#o. 



Conclusion: We build an invariant solution M on (a;^,,a;*). This solution is Lc-periodic and 
satisfies M{(jj) := {x,y, z){uj — w*), Vcj G [a;*,a;^, + Lc]- □ 
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